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(54) Title: IMAGING METHOD AND DEVICE USING SHEARING WAVES 

(54) Titre: PROCEDE ET DISPOSITIF D'IMAGERIE UTILISANT LES ONDES DE CISAILLEMENT 

(57) Abstract 

The invention concerns an imaging method 
for observing the propagation of a low-frequency 
shearing pulse wave simultaneously in multiple 
points of a diffusing viscoelastic medium (1). The 
method consists in transmitting at very high rate 
ultrasonic compression waves enabling to obtain a 
succession of images of the medium; then in delayed 
processing of the resulting images by intercorrelation 
to determine in each point of each image the 
movements of the medium while the shearing wave 
is being propagated. 

(57) Abrege- 

Proce*de* d'imagerie pour oberver la propaga- 
tion d'une onde impulsionnelle de cisaillement basse 
frequence simultan6ment en une multitude de points 
d*un milieu visco€lastique diffusant (1). A cet effet, 
on dmet a cadence ultrarapide des ondes ultrasonores 
de compression qui permettent d'obtenir une suc- 
cession d' images du milieu, puis on traite en temps 
diff6r6 les images ainsi obtenues par intercorrllation, 
pour determiner en chaque point de chaque image 
les movements du milieu lors de la propagation de „ 
l'onde de cisaillement. g „ generator 
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Precede at d i sppsit if d'imagerie utilisant ondes d^ 

cisflAllentettt- 

La presente invention est relative aux procedes et 
5 dispositifs d'imagerie utilisant les ondes de cisaille- 
ment . 

Plus particulierement , 1' invention concerne un 
procede d'imagerie utilisant les ondes de cisaillement 
pour observer un milieu viscoelastique diffusant qui con- 
10 tient des particules ref lechissant les ondes ultrasonores 
de compression, procede dans lequel on genere une onde 
elastique de cisaillement dans le milieu viscoelastique et 
on observe, au moyen d'au moins une onde ultrasonore de 
compression, le deplacement du milieu viscoelastique sou- 
15 mis a ladite onde de cisaillement. 

Le document US-A-5 810 731 decrit un exemple d'un 
tel procede, dans lequel 1 1 onde de cisaillement est gene- 
ree localement a 1 ' interieur du milieu viscoelastique ob- 
serve, au moyen de la pression de radiation d'une onde ul- 
20 trasonore modulee et focalisee sur un point a observer. On 
envoie ensuite sur ce point focal une onde ultrasonore 
supplementaire dont la reflexion permet de connaitre cer- 
tains parametres de propagation de l f onde de cisaillement 
(notamment la viscosite dynamique du milieu et son module 
25 de cisaillement) au niveau du point focal susment ionne . 

Cette technique presente 1 1 inconvenient de permet - 
tre 1' analyse d'un seul point du milieu viscoelastique 
etudie a chaque fois qu'une onde de cisaillement est gene- 
ree. Si I'on veut obtenir une image complete du milieu 
30 viscoelastique observe, il est necessaire de repeter 
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l'operation un trds grand nombre de fois, ce qui implique 
un temps de pause important (par exemple, plusieurs minu- 
tes) pour obtenir cette image. 

Ce temps de pause important rend ce proc§d£ de 
l'art anterieur peu pratique £ utiliser. 

De plus, un tel temps de pause peut nuire S 1' uti- 
lisation dudit proced£ pour obtenir une image d'un tissu 
vivant, qui est toujours en mouvement . 

La presente invention a notamment pour but de pal- 
lier ces inconvSnients . 

A cet effet, selon 1' invention, un proc£de du 
genre en question est essentiellement" caracterisi en ce 
qu'on gendre I'onde de cisaillement en appliquant au mi- 
lieu viscoelastique une excitation ayant la forme d'une 
impulsion basse frequence qui prSsente une frequence cen- 
trale f comprise entre 20 et 5000 Hz, cette impulsion 
basse frequence presentant une duree comprise entre l/2f 
et 20/f, 

en ce qu'il comporte une etape d • observation de propaga- 
tion au cours de laquelle on observe la propagation de 
I'onde de cisaillement simultanement en une multitude de 
points dans le milieu observe, ces points formant un champ 
d ■ observation sensiblement continu s'etendant au moins se- 
lon un premier axe, cette etape d • observation de propaga- 
tion de I'onde de cisaillement consistant a : 

emettre dans le milieu observe une succession 
d'au moins 10 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 
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dgtecter et enregistrer en temps reel les echos 
generes par les particules r6f lechissantes du milieu vis- 
coelastique a chaque tir d'onde ultrasonore, ces echos 
correspondant (directement ou indirectement ) a des images 
5 successives du milieu observS, 

et en ce que ledit proced§ comporte en outre une etape ul- 
tlrieure de traitement d 1 image au cours de laquelle on 
traite en temps differe les images ainsi obtenues au moins 
par intercorrelation entre images successives, pour deter- 
10 miner en chaque point du champ d' observation un param^tre 
de mouvement choisi entre le deplacement et la deformation 
du milieu visco£lastique, de fa<?on a obtenir ainsi une 
succession d' images montrant Involution du paramdtre de 
mouvement du milieu viscoelastique sous 1 ■ ef f et de la pro- 
15 pagation de I'onde de cisaillement . 

Grace a ces dispositions, on obtient un film il- 
lustrant clairement la propagation de 1 * onde de cisaille- 
ment dans le milieu viscoelastique, qui peut permettre par 
exemple, dans des applications medicales, de repgrer di- 
20 rectement des zones cancereuses dans les tissus d ! un pa- 
tient : la propagation des ondes de cisaillement s'y de- 
roule en effet tres differemment des zones voisines. 

Ce reperage s'effectue beaucoup plus facilement 
que par une observation classique par simple echographie 
25 ultrasonore, puisque la propagation des .ondes de cisaille- 
ment est fonction du module de cisaillement du milieu, 
lui-meme tres variable entre une zone de tissus sains et 
une zone de tissus cancereux : le module de cisaillement 
varie typiquement dans un rapport de 1 a 30 entre une zone 
30 saine et une zone cancereuse, alors que le module de com- 
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pression, qui regit la propagation des ondes acoustiques 
de compression utilisees dans 1 ■ echographie ultrasonore, 
varie seulement de 1 1 ordre de 5% entre un tissu sain et un 
tissu cancSreux. 

5 On notera que le film obtenu illustre la propaga- 

tion de l'onde de cisaillement beaucoup plus clairement 
que la simple succession des images donnees par les parti- 
cules ref lechissantes du milieu, puisque ledit film permet 
de visualiser £ chaque instant les zones du milieu observe 

0 qui subissent des mouvements de m§me ampleur du fait de la 
propagation de l'onde de cisaillement, alors que la suc- 
cession des images des particules ref lechissantes .ne per- 
mettrait de visualiser qu'un brouillard de points lumineux 
en mouvement . 

5 Dans des modes de realisation preferes du procSde 

selon 1' invention, on peut eventuellement avoir recours en 
outre a 1 ' une et/ou a 1 1 autre des dispositions suivantes : 
la duree de 1' impulsion basse frequence est 
comprise entre l/2f et 2f ; 

la frequence cent rale de 1' impulsion basse fre- 
quence est comprise entre 3 0 et 1000 Hz ; 

le milieu viscoelast ique observe est constitue 
par un corps vivant comprenant au moins un organe interne 
soumis a des mouvements impulsionnels , 1' impulsion basse 
frequence qui genere l'onde de cisaillement etant consti- 
tute par un mouvement impulsionnel dudit organe interne ; 

le milieu viscoelastique observe est delimite 
par une surface exterieure, et 1' impulsion basse frequence 
est appliquee au niveau de cette surface exterieure / 
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1' impulsion basse frequence est appliquee par 
un moyen d' excitation choisi parmi : 

une onde acoustique generee par au moins un trans - 
ducteur acoustique, 

et un choc genere localement par contact physique 
au niveau de la surface exterieure du milieu viscoelasti- 
que ; 

- on emet les tirs d'ondes ultrasonores de com- 
pression et on detecte les echos generis par les particu- 
les reflechissantes du milieu viscoelastique ( au moyen 
d'une batterie de transducteurs qui comprend au moins un 
transducteur et qui est disposee au contact de la .surface 
exterieure du milieu viscoelastique, 1 'onde de cisaille- 
ment etant appliquee au milieu viscoelastique en imposant 
15 un deplacement impulsionnel a ladite batterie de transduc- 
teurs ; 

ledit parametre de mouvement est la deformation 
du milieu viscoelastique : cette disposition est particu- 
lierement utile dans le dernier cas envisage ci-dessus, 
puisqu'elle permet de s'affranchir du deplacement de la 
batterie de transducteurs, deplacement qui perturberait 
sinon la mesure du deplacement des points du champ d' ob- 
servation ; 

- au cours de 1 ' observation de la propagation de 
l'onde de cisaillement , on emet entre 100 et 10 000 tirs 
d'ondes ultrasonores de compression a une cadence comprise 
entre 100 et 100 000 tirs par seconde ; 

- le champ d ' observation s • etend au moins selon 
un plan comprenant d'une part, le premier axe et d' autre 
part, un deuxieme axe perpendiculaire au premier axe ; 
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au cours de l'etape d 1 observation de propaga- 
tion, on utilise une batterie de plusieurs transducteurs 
acoustiques disposes au moins selon le deuxieme axe pour 
emettre les tirs d'ondes ultrasonores de compression et 
5 detecter les echos generes par les particules reflechis- 
santes du milieu viscoelastique , les echos d£tect6s par 
chaque transducteur acoustique etant memorises directement 
sans traitement prealable au cours de l'etape d 1 observa- 
tion de propagation, et l'etape de traitement d' image com- 
10 prenant une sous-6tape preliminaire de formation de voies 
au cours de laquelle on genere une image du milieu visco- 
elastique correspondant a chaque tir d'onde ultraspnore de 
compression, par combinaison d'au moins certains des echos 
re<?us par les diff€rents transducteurs ; 
15 - l'etape de traitement d' images est suivie (im- 

mediatement ou non) par une etape de visualisation au 
cours de laquelle on visualise au ralenti un film consti- 
tue par la succession des images traitees, chaque point de 
chaque image presentant un parametre optique qui varie 
20 suivant la valeur du parametre de mouvement affecte a ce 
point / 

le parametre optique est choisi parmi le niveau 
de gris et le niveau chromatique ; 

l'etape de traitement d' images est suivie (im- 
25 mediatement ou non) par une etape de cartographie au cours 
de laquelle, a partir de 1' evolution du parametre de mou- 
vement au cours du temps dans le champ d ' observation, on 
calcule au moins un parametre de propagation de I'onde de 
cisaillement en au moins certains points du champ d'obser- 
30 vation / 
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- le parametre de propagation de l'onde de ci- 
saillement qui est calcule au cours de l'etape de carto- 
graphie est choisi parmi la vitesse des ondes de cisaille- 
ment, le module de cisaillement , 1 ' attenuation des ondes 
de cisaillement, 1' elasticity de cisaillement, et la vis- 
cosite de cisaillement. 

Par ailleurs, 1' invention a egalement pour objet 
un dispositif d'imagerie utilisant les ondes de cisaille- 
ment pour observer un milieu viscoelastique diffusant qui 
contient des particules r£f lechissant les ondes ultrasono- 
res de compression, ce dispositif comprenant des moyens 
d' excitation pour generer une onde elastique de cisaille- 
ment dans le milieu viscoelastique et des moyens d' acqui- 
sition pour observer, au moyen d-au moins une onde ultra - 
sonore de compression, le deplacement du milieu viscoelas- 
tique soumis a ladite onde de cisaillement, 

caracterise en ce que les moyens d' excitation sont adaptes 
pour appliquer au milieu viscoelastique une excitation 
ayant la forme d'une impulsion basse frequence qui pre- 
sente une frequence centrale f comprise entre 20 et 5000 
Hz, cette impulsion basse frequence presentant une duree 
comprise entre l/2f et 20/f, 

en ce que les moyens d • acquisition sont adaptes pour ob- 
server la propagation de l'onde de cisaillement simultane- 
ment en une multitude de points dans le milieu observe, 
ces points formant un champ d 1 observation sensiblement 
continu s'etendant au moins selon un premier axe, lesdits 
moyens d ' acquisition etant adaptes pour : 

emettre dans le milieu observe une succession 
d'au moins 10 tirs d' ondes ultrasonores de compression 
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une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 

detecter et enregistrer en temps reel les echos 
generes par les particules ref lechissantes du milieu vis- 
coelastique a chaque tir d f onde ultrasonore, ces echos 
correspondant ( direct ement ou indirectement ) £ des images 
successives du milieu observe, 

et en ce que ledit dispositif comporte en outre des moyens 
de traitement d' image adaptes pour traiter en temps dif fe- 
re les images obtenues a partir des moyens d ' observation, 
au moins par intercorrelation entre images successives, 
pour determiner en chaque point du champ d' observation un 
parametre de mouvement choisi entre le deplacement et la 
deformation du milieu viscoelastique, de fa<?on a. obtenir 
ainsi une succession d'images montrant 1 1 evolution du pa- 
rametre de mouvement du milieu viscoelastique sous l'effet 
de la propagation de l'onde de cisaillement . 

Dans des modes de realisation preferes du disposi- 
tif selon 1' invention, on peut eventuellement avoir re- 
cours en outre a 1 ■ une et/ou a 1' autre des dispositions 
suivantes : 

les moyens d ' observation comprennent une batte- 
rie de transducteurs qui inclut au moins un transducteur 
et qui est adaptee pour etre disposee au contact d'une- 
surface exterieure delimitant le milieu viscoelastique, 
les moyens d 1 excitation etant adaptes pour imposer un de- 
placement impulsionnel a ladite batterie de transduc- 
teurs ; 

ledit parametre de mouvement est la deformation 
du milieu viscoelastique ; 
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- le champ d • observation s'etend au moins selon 
un plan cornprenant d'une part, le premier axe et d' autre 
part, un deuxieme axe perpendiculaire au premier axe, la 
. batterie de transducteurs comportant plusieurs transduc- 
5 teurs disposes au moins selon le deuxieme axe, des moyens 
de commande Stant prevus pour selectivement faire fonc- 
tionner le dispositif soit dans le mode d'imagerie par on- 
des de cisaillement , soit dans un mode d ■ echographie stan- 
dard permettant d'acquerir entre 10 et 100 images par se- 
10 conde . 

D'autres caracteristiques et avantages de 1' inven- 
tion apparaitront au cours de la description suiyante de 
plusieurs de ses formes de realisation, donnees a tit re 
d'exemples non limitatifs, en regard des dessins joints. 
15 Sur les dessins : 

la figure 1 est une vue schematique d'un dispo- 
sitif d'imagerie par ondes de cisaillement selon une forme 
de realisation de 1' invention, 

et la figure 2 est une vue de detail montrant 
20 une variante du dispositif de la figure 1. 

Sur les differentes figures, les memes references 
designent des elements identiques ou similaires. 

La figure 1 represente un exemple de dispositif 
d'imagerie par ondes de cisaillement selon 1 'invention, 
25 pour etudier la propagation des ondes elastiques de ci- 
saillement dans un milieu viscoelast ique 1 qui est diffu- 
sant vis a vis des ondes ultrasonores de compression, et 
qui peut etre par exemple : 
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un corps inerte, riotamment dans le cas du con- 
trole de qualite pour des applications industrielles, no- 
tamment agro- aliment aires , 

ou un corps vivant, par exemple une partie du 
corps d'un patient, dans le cas des applications medica- 
les . 

Ce dispositif comporte un transducteur acoustique 
tel qu'un haut-parleur 2 ou un pot vibreur qui est dispose 
contre la surface exterieure 3 du milieu observe 1, cette 
surface 3 Stant constitute par exemple par la peau du pa- 
tient dans les applications medicales . 

Le haut-parleur 2 peut etre command^ par un .micro- 
ordinateur 4, par exemple par 1 ' intermediaire d'un circuit 
ginerateur d f impulsions basse frequence G (ce circuit peut 
etre constitue notamment par la carte son du micro- 
ordinateur 4) et d'un amplif icateur A, pour appliquer a la 
surface 3 du milieu observe une excitation sous la forme 
d'une impulsion basse frequence, de fa<?on a generer une 
onde de cisaillement dans le milieu viscoelastique 1. 

Cette impulsion basse frequence presente en gene- 
ral une amplitude pouvant etre de 1 1 ordre de 1 mm et une 
frequence centrale f comprise entre 20 et 5000 Hz, appli- 
quee pendant une duree comprise entre l/2f et 20/f. De 
preference, la duree d 1 application de 1 'impulsion basse 
frequence est comprise entre l/2f et 2f et la frequence f 
est comprise entre 30 et 1000 Hz, cette frequence pouvant 
typiquement etre de 1 1 ordre de 50 Hz. 

En variante I'onde acoustique de cisaillement 
pourrait egalement etre obtenue (avec les caracteristiques 
susmentionnees d 1 amplitude et de frequence) : 
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par un choc genere localement par contact phy- 
sique au niveau de la surface extSrieure du milieu visco- 
elastique, par au moins un actionneur mecanique automati- 
que commande par le micro-ordinateur 4, 
5 - par un choc applique manuel lenient par contact 

physique au niveau de la surface extSrieure du milieu vis- 
coelastique, 

ou encore, par un mouvement impulsionnel natu- 
rel d'un organe interne du corps humain ou animal (par 
10 exemple un battement du coeur) dans les applications medi- 
cales. 

L'onde elastique de cisaillement produite. par le 
haut-parleur 2 se deplace avec une vitesse Cs relativement 
faible, de 1'ordre de quelques m/s (typiquement, de 1 a 10 
m/s dans le corps humain) en produisant des mouvements in- 
ternes dans le milieu viscoSlast ique 1 observe. 

Ces mouvements sont suivis en envoyant dans le mi- 
lieu 1, des ondes ultrasonores de compression qui inter- 
agissent avec les particules diffusantes 5 contenues dans 
le milieu 1, lesquelles particules sont ref lechissantes 
pour les ondes ultrasonores de compression. Les particules 
5 peuvent §tre constitutes par toute het<§rogeneite du mi- 
lieu 1, et notamment, lorsqu f il s'agit d'une application 
medicale, par des particules de collagene presentes dans 
les tissus humains . 

Pour observer la propagation de l'onde de ci- 
saillement, on utilise done une sonde ultrasonore 6 dispo- 
see contre la surface exterieure 3 du milieu observe 1. 
Cette sonde envoie, selon un axe X, des impulsions d' ondes 
ultrasonores de compression du type de celles couramment 
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utilisees en echographie, a une frequence comprise par 
exemple entre 1 et 10 0 MHz et de preference entre 3 et 15 
MHz. On notera que la sonde 6 peut etre disposee : 

soit du meme c6t£ du milieu 1 que le haut- 
5 parleur 2, comme represent^ sur la figure 1, 

soit a 1 ! oppose du haut-parleur 2 par rapport 
au milieu 1, 

soit dans toute autre position, par exemple 
dans une disposition transversale par rapport au haut- 
10 parleur 2. 

La sonde ultrasonore 6 est constitute par une bat- 
terie de n transducteurs ultrasonores Tl, T2 , . . Ti, 
. . . , Tn, n etant un nombre entier n au moins egal a 1. 

Cette sonde 6 se prtsente le plus souvent sous la 
15 forme d'une barrette linSaire pouvant comprendre par exem- 
ple n « 12 8 transducteurs alignts selon un axe Y perpendi- 
culaire £ 1 1 axe X, envoyant simultanement leurs impulsions 
d'ondes ultrasonores de facpon a generer une onde "plane" 
(c'est a dire en 1' occurrence une onde dont le front 
20 d'onde est rectiligne dans le plan X, Y) ou tout autre 
type d'onde eclairant 1' ensemble du champ d' observation . 

En variante, la batterie 2 de transducteurs peut 
eventuellement etre reduite* a un transducteur unique Tl , 
ou au contraire se presenter sous la forme d'un reseau a 
25 deux dimensions s'ttendant par exemple selon un plan per- 
pendiculaire a 1 ' axe X. 

Chacun des transducteurs Tl, T2 , ... Tn est com- 
mande par le micro-ordinateur 4 ou par une unite centrale 
CPU (qui est contenue par exemple dans une baie electroni- 
30 que 7 reliee par un cable souple a la sonde 6), pour emet- 
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tre des tirs successifs d'ondes de compression ultrasono- 
res dans le milieu 2, au cours d'une phase d ' observation 
qui peut durer par exemple moins d'une seconde et au cours 
de laquelle sont emis p tirs d'ondes ultrasonores de com- 
5 pression (p etant un entier compris entre 100 et 10 000 et 
de pr6f£rence entre 1000 et 100 000) , & une cadence com- 
prise entre 100 et 100 000 tirs par seconde et de prefe- 
rence comprise entre 1000 et 100 000 tirs par seconde, no- 
tamment entre 1000 et 10 000 tirs par seconde (cette ca- 
10 dence est limit§e par le temps d'aller-retour- de 1 ' onde de 
compression dans le milieu 1, done par l'epaisseur du mi- 
lieu 1 dans la direction X : il faut en effet que tous les 
echos generes par l'onde de compression aient ete re<?us 
par. la sonde 6 avant d'envoyer une nouvelle onde de com- 
15 pression) . 

Les tirs d'ondes ultrasonores de compression de la 
phase d' observation commencent de preference juste avant 
1' emission de 1 ' onde de cisaillement. 

De plus, dans le cas ou l'onde de cisaillement est 
20 gSneree par un mouvement impulsionnel d'un organe d'un 
corps vivant, on peut avantageusement synchroniser le de- 
marrage des tirs d'ondes ultrasonores de compression avec 
ce mouvement impulsionnel. par exemple dans le cas ou 
l'onde de cisaillement est generee par un battement de 
25 cceur, on peut synchroniser le demarrage des tirs d'ondes 
ultrasonores de compression avec une phase choisie de 
1 ' electrocardiogramme . 

Chacun de ces tirs donne lieu a la propagation 
d'une onde ultrasonore de compression dans le milieu l, 
30 avec une vitesse de propagation beaucoup plus elevee que 
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les ondes de cisaillement , par exemple de 1 ' ordre de 1500 
m/s dans le corps humain. 

L'onde ultrasonore ainsi g^neree interagit avec 
les particules r§f lechissantes 5, ce qui genere des echos 
5 ou autres perturbations analogues du signal, connus en soi 
sous le nom de "bruits de speckle" dans le domaine de 
1 1 echographie . 

Ces "bruits de speckle" sont captes par les trans - 
ducteurs TI, ... ,Tn apres chaque tir. Le signal sij (t) ain- 
10 si capte par chaque transducteur Ti apres le tir n° j est 
tout d'abord echantillonne a haute frequence (par exemple 
de 30 a 100 MHz) et numerise en temps reel (par exemple 
sur 8 bits, ou dans certains cas sur 1 bit) par un echan- 
tillonneur appartenant a la baie 7 et reli§ a ce transduc- 
15 teur, respectivement El, E2, ... En. 

Le signal sij (t) ainsi echantillonne et numerise 
est ensuite memorise, egalement en temps reel, dans une 
memoire Mi appartenant a la baie 7 et propre au transduc- 
teur Ti . 

20 Chaque memoire Mi presente par exemple une capaci- 

te de 1' ordre de 1 Mo, et contient 1 1 ensemble des signaux 
sij (t) regus successivement pour les tirs j = 1 a p. 

En temps differe, apres la memorisation de tous 
les signaux sij (t) correspondant a une meme propagation 

25 d'onde de cisaillement , l'unite centrale CPU fait retrai- 
ter ces signaux par un circuit sommateur S appartenant a 
la baie 7 (ou bien elle effectue elle-meme ce . traitement , 
ou encore ledit traitement peut etre effectue dans le mi- 
cro-ordinateur 4), par un processus classique de formation 

30 de voies . 
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On genere ainsi des signaux Sj (x, y) qui corres- 
pondent chacun a 1' image du champ d' observation apres le 
tir n° j . 

Par exemple, on peut determiner un signal Sj (t) 
5 par la formule suivante : 

SJ (0 = £ a, (x, y)Mj[t(x, y) + d,. (x, y) I V] 

OU : 

sij est le signal brut pergu par le transduc- 
teur n° i apr£s le tir d'onde ultrasonore de compression 
10 n° j, . 

t(x,y) est le temps mis par 1 ' onde ultrasonore 
de compression pour atteindre le point du champ 
d' observation de coordonnees (x,y), avec t = 0 au debut du 
tir n° j , 

15 - di(x,y) est la distance entre le point du champ 

d' observation de coordonnees (x,y) et le transducteur n° 
i, ou une approximation de cette distance, 

V est la vitesse moyenne de propagation des on- 
des acoustiques ultrasonores de compression dans le milieu 
20 viscoelastique observe. 

et ai(x / y) est un coefficient de ponderation 
tenant compte de lois d' apodisation (en pratique, on pour- 
ra dans de nombreux cas considerer que ai(x / y)=l). 

La formule ci-dessus s' applique mutatis mutandis 
25 lorsque le champ d' observation est a 3 dimensions (reseau 
plan de transducteurs) , en remplagant les coordonnees spa- 
tiales (x,y) par (x,y,z). 
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Lorsque la sonde 6 comporte un seul transducteur, 
l'etape de formation de voies est inutile, et on a direc- 
tement Sj (x) = sj[2.x/V], avec les memes notations . que 
ci-dessus . 

Apres l'etape eventuelle de formation de voies, 
1' unite centrale CPU memorise dans une mSmoire centrale M 
appartenant £ la baie 7 les signaux d' images Sj(x,y) ou 
Sj (x) ou Sj(x,y,z), qui correspondent chacun au tir n° j. 
Ces signaux peuvent egalement etre memorises dans le mi- 
cro-ordinateur 4 lorsqu'il effectue lui-meme le traitement 
d ' image . 

Ces images sont ensuite traitees deux a deux, ou- 
jours en temps differ^, par intercorrelat ion . Cette inter- 
correlation peut etre rgalisee dans un circuit DSP appar- 
tenant £ la baie 7, ou etre programmee dans 1' unite cen- 
trale CPU ou dans le micro-ordinateur 4. 

A titre d'exemple, cette intercorrelation peut se 
faire en comparant les signaux Sj (x,y) et Sj+1 (x,y) 
(dans le cas d'un champ d' observation a 2 dimensions) sur 
des fenetres spatiales glissantes de longueur Ax predeter- 
mine, pouvant aller par exemple de X a 10X, ou X est la 
longueur d'onde des ondes ultrasonores de compression 
(soit environ 0,42 a 4,2 mm a 3,5 MHz dans l'eau ou le 
corps humain) . De plus, les fenetres susmentionnees peu- 
vent se recouvrir l'une l'autre sur environ 20% de leur 
longueur selon 1 ' axe X. 

Au cours de ce processus d' intercorrelation, on 
maximise une fonction d 1 intercorrelation 

<Sj (x,y) , Sj+1 (x,y) > afin de determiner le deplacement subi 
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par chaque particule 5 dormant lieu a un echo ultrasonore, 
dans la direction X. 

Des exemples de tels calculs d'intercorrelation 
sont donnes dans l'etat de la technique, notamment par 
O'Donnell et al . (« Internal displacement and strain ima- 
ging using speckle tracking », IEEE transactions on ultra- 
sonic, f erroelectrics, and frequency control, vol. 41, n° 
3, mai 1994, p. 314-325) et par Ophir et al . 
(« Elastography : a quantitative method for imaging the 
elasticity of biological tissues », Ultrasonic imag. , vol. 
13, p. 111-134, 1991) . 

On obtient ainsi une succession de champs. de de- 
placements Djx(x,y) du milieu 1 sous l'effet de 1 ■ onde de 
cisaillement , dans la direction X. 

Cette succession de champs de deplacements est 
stockee dans la memoire M ou dans le micro-ordinateur 4 et 
peut etre visualisee, notamment au moyen de l'ecran 4a du 
micro-ordinateur, sous la forme d'un film ralenti ou la 
valeur des deplacements est illustree par un parametre op- 
tique tel que par un niveau de gris ou par un niveau chro- 
ma tique . 

On visualise ainsi parfaitement les differences de 
propagation de l'onde de cisaillement entre les zones de 
caracteristiques differentes du milieu l, par exemple les 
tissus sains et les tissus cancereux dans le cas d'une ap- 
plication medicale. 

Ce film de propagation de l'onde de cisaillement 
est en outre superposable avec une image echographique 
classique, qui peut etre generee par le dispositif decrit 
ci-dessus, lequel est apte a fonctionner : 
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soit dans le mode d'imagerie par ondes de ci- 

saillement , 

soit dans un mode d ' echographie standard, en 
fonction des commandes recues par exemple depuis le cla- 
5 vier 4b du micro-ordinateur . 

Par ailleurs, il est egalement possible de calcu- 
ler non pas les deplacements de chaque point du milieu ob- 
serve 1, mais les deformations Ejx(x,y) du milieu dans la 
direction X, c'est a dire les derivees des deplacements 
10 Djx(x,y) par rapport a x. 

Ces champs de deformations successifs sont utili- 
sables comme precedemment pour visualiser clairement la 
propagation de l'onde de cisaillement sous la forme d'un 
film, et presentent en outre 1 ' avantage de s'affranchir 
des deplacements de la sonde 6 par rapport au milieu ob- 
serve 1 . 

Cette variante est particulierement interessante 
dans la forme de realisation de la figure 2, ou la sonde 6 
est portee par le haut-parleur ou vibreur 2, ce qui impli- 
que forcement des mouvements de ladite sonde puisque c'est 
elle-meme qui genere alors l'onde de cisaillement. 

A partir des champs de deplacements ou de deforma- 
tions, on peut le cas echeant proceder en outre a une 
etape de cartographie au cours de laquelle, a partir de 
1' evolution du parametre de mouvement (deplacement ou de- 
formation) au cours du temps dans le champ d • observation 
X, Y (ou X dans la ces d'un seul transducteur , ou X, Y, Z 
dans le cas d'un reseau plan de transducteurs) , on calcule 
au moins un parametre de propagation de l'onde de ci- 
saillement, soit en certains points du champ d ' observation 
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choisis par 1 1 utilisateur a partir du micro-ordinateur 4, 
soit dans tout le champ d ' observation . 

Le paramdtre de propagation de 1 ' onde de cisaille- 
ment qui est calculi au cours de l'Stape de cartographie 
5 est choisi par exemple parmi la vitesse Cs des ondes de 
cisaillement , le module de cisaillement \i, 1 ' attenuation a 
des ondes de cisaillement, 1'Slasticite de cisaillement 
Hlviscosit<§ de cisaillement fi2 . 

Ce calcul est effectue par un processus classique 
10 d' inversion, dont un exemple est donne ci-apres dans le 
cas d'un champ d 1 observation £ deux dimensions (le meme 
processus s 1 appliquerait mutatis mutandis dans le cas d'un 
champ d* observation a une ou a trois dimensions, respec- 
tivement pour un seul transducteur Tl ou pour un reseau 
15 plan de transducteurs) . 

Dans cet exemple on fera 1 ' approximation que la 
viscosite de cisaillement \i2 est nulle et que le milieu 
est isotrope. 

L' equation d'onde donnant le vecteur deplacement D 
20 de chaque point du milieu 1 s'ecrit : 

p^ = A( M 3) (I), 
of 

ou p est la densite du milieu 1, \i est le module de ci- 
saillement (par hypothese reduit a sa partie reelle 
1' elasticity de cisaillement jil puisque la viscosite de 
25 cisaillement ji2 est supposee nulle) . 

Pour la premiere composante u du vecteur D, c'est 
a dire pour le deplacement du milieu 1 dans la direction 
X, on a done : 
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Apres transformee de Fourier temporelle et discre- 
tisation, cette equation peut etre ecrite sous la forme 
matricielle suivante, qui peut §tre ecrite pour chaque 
frequence du spectre de l'onde de cisaillement : 

B = HM (III), o\x : 

M est un vecteur de dimension (L+2 ) . (M+2 ) -4 
dont chaque composante vaut ^i m , c'est a dire la valeur 

locale du module de cisaillement en chaque point discreti- 
se de coordonnees (xi,y m ), ou 1 et m sont des entiers com- 
.pris entre 0 et respectivement L+l et M+l en eliminant les 
couples (l,m) valant (0,0), (0,M+1), (L+1,0) et (L+^M+l), 
L+2 et M+2 etant les nombres de points discretises dans 
1' image du milieu 1 respectivement selon les axes X et Y, 

B est un vecteur de dimension L.M, dont les 
composantes valent -Tn 2 .p.Ui m ou xn est la pulsation de 
l'onde de cisaillement basse frequence, p est la densite 
du milieu, Ui m est la transformee de Fourier temporelle du 
deplacement u au point de coordonnees (xi,y m ), 1 etant 
compris entre 1 et L et m etant compris entre 1 et M, 

et H est une matrice de dimension L.M lignes 
sur (L+2) . (M+2) -4 colonnes, dont les composantes sont ega- 
lement toutes connues a partir de 1' equation d'onde. 

En juxtaposant un nombre suffisant d' equations 
(III) correspondant respectivement a differentes disponi- 
bles du spectre de frequences de l'onde de cisaillement, 
on obtient une equation matricielle globale qui peut etre 
resolue par inversion matricielle, pour obtenir le vecteur 
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M, c'est a dire la valeur du module de cisaillement |i en 
tout point du champ d' observation . 

On peut alors en deduire eventuellement la valeur 
locale de la vitesse de propagation Cs de l'onde de ci- 
5 saillement en chaque point, par la formule : 




Le mode de calcul serait le meme en utilisant non 
plus les deplacements, mais les deformations du milieu ob- 
serve 1 . 
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REVENPICATIQWS 

1. ProcSde d 1 imagerie utilisant les ondes de ci- 
saillement pour observer un milieu viscoelastique diffu- 
sant (1) qui contient des particules (5) ref lechissant les 
ondes ultrasonores de compression, procede dans lequel on 
genere une onde elastique de cisaillement dans le milieu 
viscoelastique et on observe, au moyen d'au rnoins une onde 
ultrasonore de compression, le deplacement du milieu vis- 
coelastique (1) soumis a ladite onde de cisaillement, 
caracterise en ce qu'on genere 1 1 onde de cisaillement en 
appliquant au milieu viscoelastique une excitation ayant 
la forme d'une impulsion basse frequence qui presente une 
frequence centrale f comprise entre 20 et 5000 Hz, cette 
impulsion basse frequence presentant une duree comprise 
entre l/2f et 20/f, 

en ce qu'il comporte une etape d • observation de propaga- 
tion au cours de laquelle on observe la propagation de 
1 ' onde de cisaillement simultanement en une multitude de 
points dans le milieu observe, ces points formant un champ 
d ' observation sensiblement continu s'etendant au moins se- 
lon un premier axe (X), cette etape d ' observation de pro- 
pagation de l'onde de cisaillement consistant a : 

emettre dans le milieu observe une succession 
d'au moins 10 tirs d' ondes ultrasonores de compression a 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde, 

detecter et enregistrer en temps reel les echos 
generes par les particules ref lechissantes du milieu vis- 
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coelastique a chaque tir d'onde ultrasonore, ces echos 
correspondant a des images successives du milieu observe, 
et en ce que ledit precede comporte en outre une etape ul- 
terieure de traitement d' image au cours de laquelle on 
5 traite en temps dif f ere les images ainsi obtenues au moins 
par intercorrtlation entre images successives, pour deter- 
miner en chaque point du champ d' observation un parametre 
de mouvement choisi entre le deplacement et la deformation 
du milieu viscoelastique, de fa<?on a obtenir ainsi une 
10 succession d 1 images montrant 1' evolution du paramdtre de 
mouvement du milieu viscoelastique sous l'effet de la pro- 
pagation de l'onde de cisaillement . 

2. ProcedS selon la revendication 1, dans lequel 
la duree de 1' impulsion basse frequence est comprise ejatre 

15 l/2f et 2f. 

3. Procedg selon la revendication 1 ou la revendi- 
cation 2, dans lequel la frequence centrale de 1 ' impulsion 
basse frequence est comprise entre 30 et 1000 Hz. 

4. Procede selon 1 ' une quelconque des revendica- 
20 tions precedentes, dans lequel le milieu viscoelastique 

(1) observe est const itue par un corps vivant comprenant 
au moins un organe interne soumis a des mouvements impul- 
sionnels, l'impulsion basse frequence qui genere .1 ' onde de 
cisaillement etant constitute par un mouvement impulsion- 
25 nel dudit organe interne. 

5. Procgde selon 1 1 une quelconque des revendica- 
tions 1 a 3, dans lequel le milieu viscoelastique (1) ob- 
serve est delimite par une surface exterieure (3), et 
l'impulsion basse frequence est appliquee au niveau de 

30 cette surface exterieure. 
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6. Proceed selon la revendication 5, dans lequel 
1' impulsion basse frequence est appliquee par un mcyen 
d f excitation choisi parmi : 

une onde acoustique gSneree par au moins un 
transducteur acoustique (2) , 

et un choc genere localement par contact physi- 
que au niveau de la surface exterieure (3) du milieu vis- 
coelastique . 

7. Procede selon la revendication 5 ou la revendi- 
cation 6, dans lequel on emet les tirs d'ondes ultrasono- 
res de compression et on detecte les echos generes par les 
particules ref lechissantes du milieu visco^lastique, au 
moyen d'une batterie de transducteurs (6) qui comprend au 
moins un transducteur (Tl, Tn) et qui est dispos§e au 
contact de la surface exterieure du milieu viscoelast ique , 
I'onde de cisaillement etant appliquee au milieu visco- 
elastique en imposant un deplacement impulsionnel a iadite 
batterie de transducteurs. 

8. Procede selon la revendication 1, dans lequel 
ledit parametre de mouvement est la deformation du milieu 
viscoelastique . 

9. Procede selon I'une quelconque des revendica- 
tions precedentes, dans lequel, au cours de 1 1 observation 
de la propagation de I'onde de cisaillement, on emet entre 
100 et 10 000 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 1000 et 100 000 tirs par se- 
conde . 

10. Procede selon 1 ' une quelconque des revendica- 
tions precedentes, dans lequel le champ d ' observation 
s'etend au moins selon un plan comprenant d'une part, le 
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premier axe (X) et d' autre part, un deuxieme axe perpendi- 
culaire au premier axe (Y) . 

11. Precede selon la revendicat ion 10 , dans le- 
quel, au cours de l'etape d ■ observation de propagation, on 
5 utilise une batterie (6) de plusieurs transducteurs acous- 
tiques (Tl, Tn) disposes au moins selon le deuxieme axe 
(Y) pour emettre les tirs d'ondes acoustiques ultrasonores 
de compression et detecter les echos generes par les par- 
ticules r€f lechissantes (5) du milieu viscoelastique , les 
10 echos detectes par chaque transducteur acoustique etant 
memorises directement sans traitement prealable au cours 
de l'etape d ■ observation de propagation, et l'etape de 
traitement d ! image comprenant une sous-etape preliminaire 
de. formation de voies au cours de laquelle on gendre une 
15 image du milieu viscoelastique correspondant a chaque tir 
d'onde ultrasonore de compression, par combinaison d • au 
moins certains des echos reyus par les differents trans- 
ducteurs (Tl, Tn) . 

12. Procede selon l'une. quelconque des revendica- 
20 tions precedentes, dans lequel l'etape de traitement 
d' images est suivie par une etape de visualisation au 
cours de laquelle on visualise au ralenti un film consti- 
tue par la succession des images traitees, chaque point de 
chaque image presentant un parametre optique qui varie 
25 suivant la valeur du parametre de mouvement affecte a ce 
point . 

13. Procede selon la revendication 12, dans lequel 
le parametre optique est choisi parmi le niveau de gris et 
le niveau chromatique. 
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14. Procede selon l'une quelconque des revendica- 
tions precedentes, dans lequel 1'etape de traitement 
d' images est suivie par une etape de cartographie au cours 
de laquelle, a partir de Involution du parametre de mou- 
5 vement au cours du temps dans le champ d ■ observation, on 
calcule au moins un parametre de propagation de I'onde de 
cisaillement en au moins certains points du champ d' obser- 
vation. 

15. Procede selon la revendication 14, dans lequel 
10 le parametre de propagation de I'onde de cisaillement qui 

est calcule au cours de 1'etape de cartographie est choisi 
parmi la vitesse des ondes de cisaillement, le module de 
cisaillement, 1 ' attenuation des ondes de cisaillement, 
l'elasticite de cisaillement, et la viscosite de cisaille- 
15 ment . 

16. Dispositif d'imagerie utilisant les ondes de 
cisaillement pour observer un milieu viscoelastique diffu- 
sant (1) qui contient des particules (5) ref lechissant les 
ondes ultrasonores de compression, ce dispositif compre- 

20 nant des moyens d» excitation (2) pour generer une onde 
elastique de cisaillement dans le milieu viscoelastique et 
des moyens d 1 acquisition (CPU, Ti, Ei, Mi) pour observer, 
au moyen d'au moins une onde ultrasonore de compression, 
le deplacement du milieu viscoelastique soumis a ladite 

25 onde de cisaillement, 

caracterise en ce que les moyens d' excitation (2) sont 
adaptes pour appliquer au milieu viscoelastique (l) une 
excitation ayant la forme d'une impulsion basse frequence 
qui presente une frequence centrale f comprise entre 20 et 



WO 00/55616 



PCT/FR00/00599 



27 

5000 Hz, cette impulsion basse frequence presentant une 
duree comprise entre l/2f et 20/f, 

en ce que les moyens d ' acquisition (CPU, Ti, Ei , Mi) sont 
adaptes pour observer la propagation de l'onde de ci- 
saillement simultanement en une multitude de points dans 
le milieu observe, ces points formant un champ d' observa- 
tion sensiblement continu s'£tendant au moins selon un 
premier axe (X), lesdits moyens d ' acquisition etant adap- 
tes pour : 

emettre dans le milieu observe une succession 
d'au moins 10 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde, 

detecter et enregistrer en temps reel les echos 
generis par les particules rSf lSchissantes du milieu vis- 
coelastique a chaque tir d'onde ultrasonore, ces echos 
correspondant a des images success ives du milieu observe, 
et en ce que ledit dispositif comporte en outre des moyens 
de traitement d' image (CPU, S, DSP) adaptes pour traiter 
en temps differe les images dbtenues a partir des moyens 
d' observation , au moins par intercorrelation entre images 
successives, pour determiner en chaque point du champ 
d' observation un paramdtre de mouvement choisi entre le 
deplacement et la deformation du milieu viscoelastique , de 
facpon a obtenir ainsi une succession d' images montrant 
Involution du parametre de mouvement du milieu viscoelas- 
tique sous I'effet de la propagation de 1 1 onde de ci- 
saillement . 

17. Dispositif selon la revendication 16, dans le- 
quel les moyens d 1 observation comprennent une batterie (6) 
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de transducteurs (Tl, Tn) qui inclut au moins un trans - 
ducteur et qui est adaptee pour etre disposee au contact 
d'une surface exterieure (3) delimitant le milieu visco- 
elastique (l), les moyens d'excitation (2) etant adaptes 
pour imposer un deplacement impulsionnel i. ladite batterie 
(6) de transducteurs. 

18. Dispositif selon la revendication 17, dans le- 
quel ledit parametre de mouvement est la d£formation du 
milieu viscoSlastique . 

19. Dispositif selon l r une quelconque des revendi- 
cations 17 et 18, dans lequel le champ d 1 observation 
s'etend au moins selon un plan comprenant d'une part, le 
premier axe (X) et d' autre part, un deuxidme axe (Y) per- 
pendiculaire au premier axe, la batterie (6) de transduc- 
teurs comportant plusieurs transducteurs (Tl, ~., Tn) dis- 
poses au moins selon le deuxieme axe, des moyens de com- 
mande (4, CPU) Stant prevus pour selectivement faire fonc- 
tionner le dispositif soit dans le mode d'imagerie par on- 
des de cisaillement , soit dans un mode d • echographie stan- 
dard permettant d'acquerir entre 10 et 100 images par se- 
conde . 
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